A strategy towards the extraction of the Sivers function with TMD evolution 
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The QCD evolution of the unpolarized Transverse Momentum Dependent (TMD) distribution 
functions and of the Sivers functions have been discussed in recent papers. Following such results we 
reconsider previous extractions of the Sivers functions from semi-inclusive deep inelastic scattering 
data and propose a simple strategy which allows to take into account the Q 2 dependence of the 
TMDs in comparison with experimental findings. A clear evidence of the phenomenological success 
of the TMD evolution equations is given, mostly, by the newest COMPASS data off a transversely 
polarized proton target. 
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I. INTRODUCTION AND FORMALISM 



The exploration of the 3-dimensional structure of the nucleons, both in momentum and configuration space, 
is one of the major issues in hadron high energy physics, with dedicated experimental and theoretical efforts. 
In particular, several Semi-Inclusive Deep Inelastic Scattering (SIDIS) experiments are either running or being 
planned. From the measurements of azimuthal asymmetries, both with unpolarized and polarized nucleons, one 
obtains information on the Transverse Momentum Dependent Parton Distribution Functions (TMD PDFs) and 
on the Transverse Momentum Dependent Fragmentation Functions (TMD FFs). The TMD PDFs and the TMD 
FFs are often globally referred to simply as TMDs. The TMD PDFs convey information on the momentum 
distributions of partons inside protons and neutrons. 

The analysis of the experimental data is based on the so-called TMD factorization, which links measurable 
cross sections and spin asymmetries to a convolution of TMDs. In particular, the Sivers function, which 
describes the number density of unpolarized quarks inside a transversely polarized proton, has received much 
attention and has been extracted from SIDIS data by several groups, with consistent results [TH5]. However, 
all these phenomenological fits of the Sivers function (and other TMDs) have been performed so far using a 
simplified version of the TMD factorization scheme, in which the QCD scale dependence of the TMDs - which 
was unknown - is either neglected or limited to the collinear part of the unpolarized PDFs. While this might not 
be a serious numerical problem when considering only experimental data which cover limited ranges of low Q 2 
values, it is not correct in principle, and taking into account the appropriate Q 2 evolution might be numerically 
relevant for predictions at higher Q 2 values, like future electron-ion or electron-nucleon colliders (EIC/ENC) 
and Drell-Yan experiments. 

Recently, the issue of the QCD evolution of unpolarized TMDs and of the Sivers function has been studied in a 
series of papers [7H9] and a complete TMD factorization framework is now available for a consistent treatment of 
SIDIS data and the extraction of TMDs. A first application of the new TMD evolution equations to some limited 
samples of the HERMES and COMPASS data [TO] has indeed shown clear signs of the Q 2 TMD evolution. 

We follow here Refs. [5] and [5] adopting their formalism, which includes the explicit Q 2 dependence of the 
TMDs, and apply it to the extraction of the Sivers function from SIDIS data, exploiting the latest HERMES [11 



and COMPASS [Tj?] results. In the sequel of this Section we present the explicit formalism: in Subsection I A 



we describe the setup and structure of the TMD evolution equations, in Subsection |I B| we discuss the parame 
terizations used for the unknown input functions, while in Subsection I C we present analytical solutions of the 
TMD evolutions equations obtained under a specific approximation. 

In Section |h] we perform a best fit of the SIDIS Sivers asymmetries taking into account the different Q 2 values 
of each data point and the Q 2 dependence of the TMDs; we compare our results with a similar analysis performed 
without the TMD evolution. Differences between Sivers functions extracted from data with and without the 
TMD evolution are shown and commented. In all this we differ from Ref. [10], which explicitly shows the 
evolution of an existing fit of the Sivers SIDIS asymmetry [T3] from the average value (Q 2 ) = 2.4 GeV 2 for 
HERMES data [HJ to the average value of (Q 2 ) = 3.8 GeV 2 for the most recent COMPASS data [JJ. Further 
comments and conclusions are given in Section [HI| 
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A. Formalism for TMD Q 2 dependence 

In Refs. [7] and [5], Collins, Aybat and Rogers have proposed a scheme to describe the Q 2 evolution of the 
TMD unpolarized distribution and fragmentation functions: within the framework of the Collins-Soper-Sterman 
(CSS) factorization formalism |14|I15]. they can describe the non-perturbative, low transverse momentum region 
and, at the same time, consistently include the perturbative corrections affecting the region of larger energies 
and momentum transfers. However, this formalism cannot be directly applied to spin dependent distribution 
functions, like the Sivers function [T^], for which the collinear limit does not exist. 

More recently, an extension of the unpolarized TMD-evolution formalism was presented in Ref. j9] to provide 
a framework in which also spin-correlated PDFs can be accounted for. For our purposes, we will use Eq. (44) of 
Ref. [9] which, compared to the unpolarized TMD evolution scheme, Eq. (26) of Ref. jS], requires the extra aid of 
a phenomenological input function embedding the missing information on the evolved function, that, in the case 
of the Sivers function, is both of perturbative and non-pertubative nature. Although the unpolarized PDF and 
FF TMD evolution equations are in principle known [8], in this paper we adopt the simplified functional form 
of the evolution equation, as proposed for the Sivers function in Ref. [5], for all TMD functions, for consistency. 

Thus, we strictly follow Ref. [9] and combine their Eqs. (44), (43) and (30), taking, as suggested [9], the 
renormalization scale fi 2 and the regulating parameters Cf and £d all equal to Q 2 . Then, the QCD evolution 
of the TMDs in the coordinate space can be written as 



F(x, b T ; Q) = F{x, b T ; Q ) exp { In K{b T ; Q 




(1) 



where F can be either the unpolarized parton distribution, F(x, bx; Q) = f q / p (x, bx; Q), the unpolarized frag- 
mentation function F(x,bx',Q) = Dh/q{z,bT;Q), or the first derivative, with respect to the parton impact 
parameter b^, of the Sivers function, F(x,bT',Q) = /(^ (x, 6t;Q). Notice that throughout the paper in- 



dependent distribution and fragmentation functions will be denoted with a 
In the above equation the function K is given in general by [9] : 



on top. 



K{br,n) = Mb) + 



dfj, 



9k (br) , 



(2) 



with, at 0(a s ) Pi], 



K(h,fi b ) = - [\n(bl 4) - In 4 + 2 lE ] 



(3) 



K(b T ) EE - hT ^ = ^ • (4) 

The first two terms in Eq. ^ are perturbative and depend on the scale /i through the coupling a s (/x), while 
the last term is non-perturbative, but scale independent. C\ is a constant parameter which can be fixed to 
optimize the perturbative expansion, as explained in Ref. [15]. Refs. [8] and [9] adopt the particular choice 
C\ = 2e~ lE which automatically implies K(b*,/ib) = 0, considerably simplifying the bx dependence of the CSS 
kernel K{b Tl fJ.), Eq. jij. 

The anomalous dimensions and 7^ appearing respectively in Eqs. and ([2]), are given, again at order 
0(a s ),by|tt3] 

Q\ , , C F /3 . Q 2 \ , , f , 2C F 

7 ) = flsM 7ir ln ? 



7f(m; ty ) = OL a (ji) - - In — ) Jk(h) = ol s (jj) ■ (5) 



By making use of Eqs. ([2j)-([5|), the evolution of F(x, bj-', Q) in Eq. (Ilj) can then be written as 



F(x, b T ; Q) = F(x, b T - Q ) R(Q, Q , b T ) exp \ ~g K (b T ) In ^ } , (6) 
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with 



i?(Q,Q ,6 T )^exp{ln-^ f^M+ f — If ( /*, ^ )) • (7) 



The Q 2 evolution is driven by the functions gxibr) and R(Q,Qo,br)- While the latter, Eq. ([7j), can be easily 
evaluated, numerically or even analytically, the former, is essentially unknown and will need to be taken from 
independent experimental inputs. 

The explicit expression of the TMDs in the momentum space, with the QCD Q 2 dependence, can be obtained 
by Fourier-transforming Eq. ([6]), obtaining [9]: 

fq/p(x,kx;Q) = — / db T b T J (kj_b T ) fq/ p (x,b T ;Q) (8) 
Z7T Jo 

Dh/q(*>P±;Q) = ir [ db rb T J (k T b T ) D h/q (z,b T ;Q) (9) 
Z7T Jo 

-» —1 f°° ~ 

ftf{x, kr, Q) = ~z r~ / db T b T Jt(k ± b T ) f£*(x, b T] Q) , (10) 

where Jo and J\ are Bessel functions. In this paper we denote the distribution and fragmentation functions 
which depend on the transverse momenta (TMDs) with a "widehat" on top. f q / p is the unpolarized TMD 
distribution function for a parton of flavor q inside a proton, and Dh/ q is the unpolarized TMD fragmentation 

function for hadron h inside a parton q. is the Sivers distribution defined, for unpolarized partons inside a 
transversely polarized proton, as: 

f q/p t(x,k±,S;Q) = fq/ p (x,kx;Q) - At fc -U Q)^^ — (11) 

= f q/p (x,k ± ;Q) + -A N f q/p ,(x,k ± ;Q) 11 ^ ■ (12) 

In our notation fcj^ is the transverse momentum of the parton with respect to the parent nucleon direction 
and p_|_ is the transverse momentum of the final hadron with respect to the parent parton direction. Notice 
that in Refs. [5] and [S] all transverse momenta are defined in a unique frame, the so-called hadron frame, in 
which the measured hadrons have zero transverse momentum. In this frame, the initial and the final parton 
transverse momenta are denoted, respectively, by kir and k2T- They are related to our notation by: k± = ki^ 
and, at leading order in p±, = ~z\lit. This requires some attention when dealing with the fragmentation 
functions. Usually, the TMD FFs are defined in terms of the hadronic p±, i.e. the transverse momentum of the 
final hadron h with respect to the direction of the fragmenting parton q, while, following Refs. [S] and [5], the 
Fourier transform Q is performed from the impact parameter space of the fragmenting parton (bx) into the 
corresponding partonic transverse moment um ( kj = pj_/z) in the hadron frame. This will generate some extra 
z 2 factors, as explained in detail in Section I B 



B. Parameterization of unknown functions 



Eqs. QSh— QIOJ) can be adopted as the appropriate functional forms, with the correct Q 2 dependence induced 
by Eqs. ([6])-([7])j to be used in the extraction of phenomenological information on the unpolarized and Sivers 
TMDs. In order to do so, one should start with a parameterization of the unknown functions inside Eq. (jfjj): 
gnibr) and F(x,bT',Qo)- As already anticipated, gxibr) is a non-perturbative, but universal function, which 
in the literature is usually parameterized in a quadratic form. As in Refs. [3] and |10j . we will adopt the results 
provided by a recent fit of Drell-Yan data [I7j , and assume 

gxibr) = x gib\ with g 2 — 0.68 corresponding to 6 max = 0.5 GeV -1 . (13) 
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We should now parameterize the function F(a;, 6 T ;Q ) m configuration space. We wish to test the effect 
of the TMD evolution in the extraction of the Sivers functions from data; in particular we will compare the 
extraction based on TMD evolution with previous extractions which did not take such an evolution into account. 
Then, we parameterize the input function F(x, bp] Qo) by requiring that its Fourier-transform, which gives the 
corresponding TMD function in the transverse momentum space, coincides with the previously adopted k±- 
Gaussian form, with the x dependence factorized out. That was also done in Refs. [5] and [5], assuming for the 
unpolarized TMD PDF 

fq/ P (x,b T ;Q Q ) = f q/p (x,Q Q )exp{-a 2 b? r } , (14) 

where f q / p (x, Q ) is the usual integrated PDF of parton q inside proton p, evaluated at Q ; the value of a 2 is 
fixed by requiring the desired behavior of the distribution function in the transverse momentum space at the 
initial scale Qo: taking a 2 — (fcj_)/4 one recovers 

f q / P (x, k±; Qq) = f q/p (x, Qo) e ~ k ^^ , (15) 

in agreement with Refs. [5J [T21 [T5] . 

Similar relations hold for the TMD FFs, with an additional z 2 factor due to the fact that the Fourier- 
transform ([9| leads from the impact parameter space of the fragmenting parton in the hadron frame to the 
corresponding partonic transverse momentum kr, while the TMD FFs are functions of the transverse momentum 
pj_ = zkx of the final hadron with respect to the fragmenting parton direction. This requires the initial 
parameterization 

D h / q {z, b T ; Qo) = \ D h/q (z, Q ) exp {~/3 2 b^} , (16) 

where D h / q (z,Qo) is the usual integrated FF evaluated at the initial scale Qo, and /3 2 — (p^/Az 2 in order to 
recover the previously adopted behavior [51 ll3lHB] of the fragmentation function in the p± transverse momentum 
space at Q : 

D h/q (z, P± ;Qo) = D h/q (z,Q ) : ^e-rt/^ . (17) 

Analogously, we parameterize the Sivers function at the initial scale Qo as 

f&(x, b T ; Qo) = -2 7 2 ft T {x; Qo) b T e'^ h - , (18) 



which, when Fourier-transformed according to Eq. (10), yields: 



ft T {x, k x ;Q ) = ft T (x- Qo) — ^ e - fc W . (19) 

4 7T 7 



Eq. (19) agrees with our previous parameterization of the Sivers function, at the initial scale Qo [5l 1131 fT8] . 
taking: 

47 2 ^(fci)5= (20) 



firfaQo) = -^[-/ 2e ^f q/pt (x,Q ) ■ 



(21) 



Mi is a mass parameter, M p the proton mass and A"L p t (x, Qo) is the x-dependent term of the Sivers function, 
evaluated at the initial scale Qo and written as [13l [T8] : 

A N f q/pf (x,Q ) = 2N q {x) f q/p (x,Q ) , (22) 

where Af q (x) is a function of x, properly parameterized (we will come back to details of the Sivers function 
parameterization in Section \Uh. 
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The final evolution equations of the unpolarized TMD PDFs and TMD FFs, in the configuration space, are 
obtained inserting Eqs. (14) and (161 into Eq. ^6|: 



fq/p(x,b T ;Q) = f g /p(x,Q ) R(Q,Q a ,b T ) exp \ -b~ T \ n 



2 Q 



D h / q (z,b T ;Q) = -^D h/q (z,Q ) R(Q,Q ,b T ) exp< -b T [0 



,2 , 92 , Q 

T Qo 



(23) 
(24) 



with a 2 = (£rp/4, ,3 2 = (p 2 L )/(4z 2 ), 52 given in Eq. (13] and R{Q,Q ,b T ) in Eq. ([7j. 

The evolution of the Sivers function is obtained through its first derivative, inserting Eq. ( 18 1 into Eq. ([6| 



(\ ?lr. h r . ()) = -2 - - J u (,■: ()„) fn<). (),,. !„)!>, oxp 1 ' -l>j- ( 7 2 + f In ^ 



(25) 



with 7 and fpr(x; Qg) given in Eqs. ( 20 )- ( 22 ) 



Eqs. (23)-(25) show that the Q 2 evolution is controlled by the logarithmic Q dependence of the bx Gaussian 
width, together with the factor R(Q, Qo, bx)- for increasing values of Q 2 , they are responsible for the typical 
broadening effect already observed in Refs. [H] and [5]. 

It is important to stress that although the structure of Eq. ([T]) is general and holds over the whole range 
of bx values, the input function F(x,bx,Qo) is only designed to work in the large-6r region, corresponding 
to low k± values. Therefore, this formalism is perfectly suitable for phenomenological applications in the 
kinematical region we are interested in, but the parameterization of the input function should be revised in the 
case one wishes to apply it to a wider range of transverse momenta, like higher Q 2 processes where perturbative 
corrections become important. 



C. Analytical solution of the TMD Evolution Equations 



The TMD evolution in Eqs. ( 23 )-( 25 1 implies, apart from the explicit Gaussian dependence, a further non 
trivial dependence on the parton impact parameter bx through the evolution kernel R(Q, Qo, bx) and the upper 
integration limit /if,, Eq. Q, which appears in Eq. ([7]); consequently it needs to be evaluated numerically. 
However, the evolution equations can be solved analytically by making a simple approximation on this bx 
dependence. A close examination of Eq. Q shows that /i& is a decreasing function of bx that very rapidly freezes 
to the constant value Ci/6 max = /ib{bx ~> oo): more precisely, the approximation ^ = const, holds for any 
bx > 1 GeV -1 . As very small values of bx correspond to very large values of this approximation is safe in our 

framework, where the typical k± are less than 1 GeV. Neglecting the &T-dependence of the factor R(Q, Qo, bx) 
does not depend on bx anymore, see Eq. Q, and can even be integrated analytically by using an explicit 
representation of a s (Q). In the sequel we will refer to it as R(Q,Qo), with R{Q,Qo) = R{Q,Qo,bx — > oo)- 
Fig. [ljshows the evolution factor R(Q, Qo, bx) plotted as a function of bx at two fixed values of Q 2 (left panel), 
and R(Q, Qo) as a function of Q 2 (right panel). It is clear that R{Q, Qo) settles to a constant value for bx > 1 
GeV" 1 . In both cases, Q 2 , = 1 GeV 2 . 

Thus, in this approximation, the TMD evolution equation ^ only depends on bx through the non 



perturbative function gx(bx), which has been chosen to be a quadratic function of bx, Eq. (13), and through 
the bx dependence of the initial input function F(x,bx;Qo) which has been chosen to be Gaussian. It results 
in a fo^-Gaussian form, with a width which depends logarithmically on Q/Qo, for the TMD evolution equation. 
For the unpolarized TMD PDFs one has 

f q/p (x,b T ;Q) = f q/p (x,Q a ) R(Q,Q ) cxp j-y ) + 2 ff2 In ^) j . (26) 

Its Fourier-transform, Eqs. ([8]), delivers a Gaussian distribution in the transverse momentum space as well: 

e -k 2 ± /w 2 

f q /p(x, k ± ;Q) = f q/p (x, Qo) R(Q, Qo) 3— , (27) 
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FIG. 1: In the left panel, the evolution factor R(Q, Qo, br) is plotted as a function of br at two fixed values of Q 
the right panel we show R(Q, Qo) = R(Q, Qo, &t — > oo) as a function of Q 2 . In both cases, Qo = 1 GeV 2 . 



In 



where f q / p {x,Qo) is the usual integrated PDF evaluated at the initial scale Qo and, most importantly, 
w 2 (Q,Qq) is the "evolving" Gaussian width, defined as: 



w 2 (Q,Q ) = (ki)+2g 2 \n 



Qo 



(28) 



It is worth noticing that the Q 2 evolution of the TMD PDFs is now determined by the overall factor R(Q, Qo) 
and, most crucially, by the Q 2 dependent Gaussian width w(Q, Qo). 
The TMD FFs evolve in a similar way, Eq. (pIL 



D h / q (z,b T ;Q) = -^D h / q (z,Qo) R(Q,Qo) exp 



leading to the TMD FF in momentum space, 



4^2 



D h/q (z,p±;Q) = D h/q (z,Q ) R(Q,Q ) 



(pi) +2z 2 g 2 \n 



(29) 



(30) 



with an evolving and z-dependent Gaussian width wf = Wf(Q, Qo) given by 

w 2 = w 2 {Q, Qo) - (pi) + 2z 2 g 2 In % ■ 

Wo 



(31) 



For the Sivers distribution function, by Fourier-transforming Eq. ( 25 ) (with R — > R) as prescribed by Eq. ( 10 ) , 
we obtain [see also Eqs. (fTTl), Q, Q and pTJ]: 



A N f q/p r(x,k ± ;Q) 



f2e ^1 A N f q / p t (x, Qo) R(Q, Qo) 



-kl/w% 



(32) 



with 



w 2 s (Q,Qo) = (fci) s + 2 52 ln 



(33) 



It is interesting to notice that the evolution factor R(Q,Qo), controlling the TMD evolution according to 
Eqs. ([27]), (|30| and ([32]) is the same for all functions (TMD PDFs, TMD FFs and Sivers ) and is flavor 
independent: consequently it will appear, squared, in both numerator and denominator of the Sivers azimuthal 
asymmetry and, approximately, cancel out. Therefore, we can safely conclude that most of the TMD evolution 
of azimuthal asymmetries is controlled by the logarithmic Q dependence of the k± Gaussian widths w 2 (Q, Qo), 
Eqs. ([28]), ([31]) and ([33]). We will come back to this in Section [n] 



7 




k ± [GeV] k ± [GeV] 



FIG. 2: The left panel shows the unpolarized TMD PDF, f u/p , evolved from the initial scale, Q 2 , = 1 GeV 2 , to Q 2 = 2.4 
GeV 2 , using TMD-evolution (red, solid line), DGLAP-evolution (blue, dashed line) and the analytical approximated 
TMD-evolution (green dot-dashed line). The right panel shows the same functions at the scale Q 2 = 20 GeV 2 . Notice 
that, while there is hardly any difference between the DGLAP-evolved lines at Q 2 = 2.4 and Q 2 = 20 GeV 2 , the TMD 
evolution induces a fast decrease in size of the TMD PDF functions at large Q 2 and a simultaneous widening of its 
Gaussian width. Here the analytical approximated evolution gives results in good agreement with the exact calculation 
even at large Q 2 . 




FIG. 3: The left panel shows the ratio Sivers/PDF, A f u / p f/2f u / p , evolved from the initial scale, Qq = 1 GeV 2 , 
to Q 2 = 2.4 GeV 2 , using TMD-evolution (red, solid line), DGLAP-evolution (blue, dashed line) and the analytical 
approximated TMD-evolution (green dot-dashed line). The right panel shows the same functions at the scale Q 2 = 20 
GeV 2 . Notice that, while there is almost no difference between the DGLAP-evolved lines at Q 2 = 2.4 and Q 2 — 20 GeV 2 , 
the TMD evolution induces a fast decrease in size of the ratio Sivers/PDF functions with growing Q 2 and a simultaneous 
widening of its Gaussian width. It is interesting to point out that the analytical approximation, for the Sivers function, 
visibly breaks down at large values of Q 2 . 



To illustrate the features of this new TMD evolution, we compare it with the results obtained evolving only 
the collinear part, f q / p (x,Q), of the unpolarized TMD PDF according to the usual DGLAP equations and 
assuming the k± dependent term of this function to be unaffected by evolution. In the left panel of Fig. [2] we 
show the k± behavior of the unpolarized TMD PDF f u /p(x,k±,Q 2 ), at the fixed value x — 0.1, evaluated at 
the scale Q 2 = 2.4 GeV 2 (the average Q 2 value for the HERMES experiment). In the right panel we show 
the same function at a higher scale, Q 2 — 20 GeV 2 (which is the highest bin average Q 2 detected in the 
COMPASS experiment). In both cases the chosen initial scale is Q 2 , — 1 GeV 2 . The red, solid line corresponds 
to the k± distribution of the TMD PDF found by using the TMD-evolution of Eq. ( 23 ) while the blue, dashed 
line represents the result obtained by using DGLAP evolution equations. At the initial scale, Qq = 1 GeV 2 , 
solid and dashed curves coincide, by definition. However, while the DGLAP evolution is so slow that there 
is hardly any difference between the DGLAP-evolved lines at Q 2 = 2.4 GeV 2 and Q 2 = 20 GeV 2 , the TMD 
evolution induces a fast decrease of the maximum values of the TMD PDF function with growing Q 2 , and a 
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simultaneous broadening of its Gaussian width, as observed in Refs. [H) and [5]. It is interesting to notice that 
the approximated evolution of Eq. ( |27| ) , corresponding to the green, dot-dashed line works really well, even for 
large Q 2 values. 

A similar study is performed in Fig. [3] for the Sivers function. Here, by DGLAP evolution we mean that 
the Sivers function evolves like an unpolarized collinear PDF, only through the factor f q / p (x,Q) contained 
in its parameterization, Eq. (22 1. The parameters used for the plots are those given in Table [Tlj although 
any set of realistic parameters would lead to the same conclusions. The left panel shows the ratio between 
the Sivers function and the TMD PDF, A N f u / p t(x,k±;Q)/(2f u / p (x,kj_;Q)), evaluated at the scale Q 2 = 2.4 
GeV 2 . Again, the red, solid line is obtained using the TMD-evolution of Eqs. (23) and ( f25| , while the blue, 
dashed line is given by the DGLAP-evolution. The green dot-dashed line represents the results obtained using 
the approximated analytical TMD-evolution of Eqs. (27) and (32). The right panel shows the same functions 
at the scale Q 2 = 20 GeV 2 . Similarly to the case of TMD PDFs, while there is no difference between the 
DGLAP-evolved lines at Q 2 = 2.4 and Q 2 = 20 GeV 2 , the TMD evolution induces a fast decrease in the size of 
the TMD Sivers functions with growing Q 2 and a simultaneous widening of its Gaussian width. It is interesting 
to point out that the analytical TMD approximation, for the Sivers function visibly breaks down for large values 
of k±. 



II. SIDIS DATA AND TMD VS. NON-TMD EVOLUTION 



Having established the phenomenological formalism necessary to implement the TMD evolution, as given in 

N ? ■ -(x,k x ,Q) = (-2k ± /M p )ft T 



Refs. [THHj, we apply it to the Sivers function. This TMD distribution, A N f q / p t 



can be extracted from HERMES and COMPASS £p 
defined as 



.4 



sin(<p h -c/>s) 
UT 



h X SIDIS data on the azimuthal moment A 



[da^ - da 4 -] sm(<f> h - 4>s) 



J dcf) S d<t> h [da^ + cM] 



I IT' 

sin(</>h-0 s ) 
UT ' 



(34) 



This transverse single spin asymmetry (SSA) embeds the azimuthal modulation triggered by the correlation 
between the nucleon spin and the quark intrinsic transverse momentum. The "weighting" factor sin(0^ — <fis) 
in Eq. (34) is appropriately chosen to single out, among the various azimuthal dependent terms appearing in 
[da^ — da^\ , only the contribution of the Sivers mechanism |181 I19j . By properly taking into account all intrinsic 
motions this transverse single spin asymmetry can be written as pQ 



A 



l(4>h-<t>s) 



^ J d(j) S d(j) h d 2 k 1 _ A N f q/pt (x, k_i_,Q) sin(tp - <f> s ] 



da^ e i 



dQ 2 



D^(z,p±,Q) sh 



UT 



^2 / d(j) S d(j} h d 2 k_ L f q / p (x, kj_, Q) 



dQ 2 



D*(z,P±,Q) 



(35) 



With respect to the leptonic plane, 4>s and 4>h are the azimuthal angles identifying the transverse directions of 
the proton spin S and of the outgoing hadron h respectively, while ip defines the direction of the incoming (and 
outgoing) quark transverse momentum, k± = fcj_(cos ip, sin if, 0); da ILq ~^ tq j dQ 2 is the unpolarized cross section 
for the elementary scattering tq — > Iq. 

The aim of our paper is to analyze the available polarized SIDIS data from the HERMES and COMPASS 
collaborations in order to understand whether or not they show signs of the TMD evolution proposed in Ref. [S] 
and described in Section I A Our general strategy is that of adopting the TMD evolution in the extraction of 
the Sivers functions, with the same parameterization and input functions as in Refs. [5"}I13|. and see if that can 
improve the quality of the fits. In doing so we will make use of the HERMES re-analysis of SIDIS experimental 
data on Sivers asymmetries for pion and kaon production and the newest SIDIS COMPASS data off a proton 
target, which cover a wider range of Q 2 values, thus giving a better opportunity to check the TMD evolution. 

In particular we perform three different data fits: 



• a fit (TMD-fit) in which we adopt the TMD evolution equation discussed in the Section I A and IB 
Eqs. ([23j)-([25j) and ([8|-(|l0]); 



• a second fit (TMD-analytical-fit) in which we apply the same TMD evolution, but using the analytical 
approximation discussed in Section II CI Eqs. (p7|), pOl and ((32|); 
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• a fit (DGLAP-fit) in which we follow our previous work, as done so far in Ref. [SJ [T3], using the DGLAP 
evolution equation only in the collinear part of the TMDs. 

As a result of the fit we will have explicit expressions of all the Sivers functions and their parameters. However, 
the goal of the paper is not that of obtaining a new extraction of the Sivers distributions, although in the sequel 
we will show, for comment and illustration purposes, the Sivers functions for u and d valence quarks, with the 
relative parameters. The procedure followed here aims at testing the effect of the TMD evolution, as compared 
with the simple DGLAP evolution so far adopted, in fitting the TMD SIDIS data. If it turns out, as it will, 
that this improves the quality of the fit, then a new extraction of the Sivers distributions, entirely guided by 
the TMD evolution, will be necessary. That will require a different approach from the very beginning, with 
different input functions and parameterizations. 

Here, we parameterize the Sivers function at the initial scale Qq — 1 GeV, as in Ref. [SJ [T3] , in the following 
form: 



with 



A f q /pt(x, k ± ,Q ) = A f q/p t(x,Q )h(k ± ) = 2Af q (x)h(k ± ) f q/p (x,k ± ,Q ), (36) 



U q (x) = N q x aq (1 - xf* [aq+ : q [ , (37) 

h{k±) = V2e%-e- k V M ? , (38) 



Mi 



where f q / p (x, k±,Q ) is defined in Eq. 15 and N q , a q , 8 q and M% (GeV) are (scale independent) free parameters 
to be determined by fitting the experimental data. Since h(k±) < 1 for any k± and |A/" q (x)| < 1 for any x 
(notice that we allow the constant parameter N q to vary only inside the range [—1,1]), the positivity bound for 
the Sivers function, 

|A~W(x,*J, £ii 



is automatically fulfilled. Similarly to PDFs, the FFs at the initial scale are parameterized with a Gaussian 



shape, Eq. (17 1. 

As in Refs. [20] and [5] , the average values of k± and p± are fixed as 

(k 2 ± ) = 0.25 GeV 2 (pi) = 0.20 GeV 2 . (40) 

We take the unpolarized distributions f q / p (x, Qq) from Ref. |21| and the unpolarized fragmentation functions 
Dh/ q ( z ,Qo) n ' om Ref. [55], with Qq = 1.0 GeV. As in Ref. [5], we adopt 11 free parameters: 



N Uv 


N dl 


N s 






N-s 




a d v 







Mi 


(GeV) , 



(41) 

where the subscript v denotes valence contributions. In this choice we differ from Ref. [5], where valence and 
sea contributions were not separated. 

We perform best fits of 11 experimental data sets: HERMES [11 data for SIDIS production of pions (tt + , 
7r~, 7T°) and kaons (K + and K~), COMPASS data for SIDIS pion (7r + , tt~) and kaon (K + and K~) production 
from a LiD (deuteron) target [53], and the preliminary COMPASS data for charged hadron production from 
an NH3 (proton) target [T5J. The results of these 3 fits are presented in Table [i] in terms of their x 2 s. 

As it is clear from the first line of Table I, the best total Xtot> which amounts to 256, is obtained by using 
the TMD evolution, followed by a slightly higher xiot °f the analytical approximation, and a definitely larger 
Xtot — 316 corresponding to the DGLAP fit. To examine the origin of this difference between TMD and DGLAP 
evolution, we show the individual contributions to xiot °f eacn experiment (HERMES, COMPASS on NH 3 and 
on LiD targets), for all types of detected hadrons and for all variables observed (x, z and Pt)- A global look 
at the numbers reported in Tabic I shows that the difference of about 60 x 2 -points between the TMD and the 
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TABLE I: x 2 contributions corresponding to the TMD-fit, the TMD-analytical-fit and the DGLAP-fit, for each experi- 
mental data set of HERMES and COMPASS experiments. 
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FIG. 4: The results obtained from our fit of the SIDIS J 4^y''* h ~'* s ' Sivers asymmetries applying TMD evolution (red, 
solid lines) are compared with the analogous results found by using DGLAP evolution equations (blue, dashed lines). 
The green, dash-dotted lines correspond to the results obtained by using the approximated analytical TMD evolution 
(see text for further details). The experimental data are from HERMES [TT] (left panel) and COMPASS [H] (right 
panel) Collaborations. 



DGLAP fits is not equally distributed among all x 2 s per data point; rather, it is heavily concentrated in three 
particular cases, namely in the asymmetry for tt + production at HERMES and for h + and h~ production at 
COMPASS off a proton target, especially when this asymmetry is observed as a function of the x- variable. 

It is important to stress that, as x is directly proportional to Q 2 through the kinematical relation Q 2 = 
xy s, the x behavior of the asymmetries is intimately connected to their Q 2 evolution. While the HERMES 
experimental bins cover a very modest range of Q 2 values, from 1.3 GeV 2 to 6.2 GeV 2 , COMPASS data raise 
to a maximum Q 2 of 20.5 GeV 2 , enabling to test more severely the TMD Q 2 evolution in SIDIS. 

These aspects are illustrated in Fig. [4J where the SIDIS Sivers asymmetries A s ^:^ h ~^ s " > obtained in the three 
fits are shown in the same plot. It is evident that the DGLAP evolution seems to be unable to describe the 
correct x trend, i.e. the right Q 2 behavior, while the TMD evolution (red solid line) follows much better the 
large Q 2 data points, corresponding to the last ir-bins measured by COMPASS. The approximate analytical 
TMD evolution (green dash-dotted line) works very well for low to moderate values of Q 2 while it starts to 
deviate from the exact behavior at large Q 2 values. 

In Figs. [5] we show, as an illustration of their qualities, our best fits (solid red lines) of the HERMES 
experimental data IV on the Sivers asymmetries for pion production. Those on the left panels are obtained 
adopting the new TMD evolution, while those on the right use the simplified DGLAP evolution. Similar results 
are shown, for the recent COMPASS data off a proton target [12] for charged hadron production, in Fig. [6] 

The shaded area represents the statistical uncertainty of the fit parameters corresponding to a Ax 2 — 20 (i.e. 
to 95.45% confidence level for 11 degrees of freedom, see Appendix A of Ref. [5] for further details). Notice 
that, in general, the error bands corresponding to the TMD-evolution fit are thinner than those corresponding 
to the DGLAP fit: this is caused by the fact that the TMD evolution implies a ratio Sivers/PDF which becomes 
smaller with growing Q 2 , as shown in Fig. [3J constraining the free parameters much more tightly than in the 
DGLAP-evolution fit, where the Sivers/PDF ratio remains roughly constant as Q 2 raises from low to large 
values. 

In Fig. [7] we compare, for illustration purposes, the Sivers function - actually, its first moment, defined in 
Ref. [5 - at the initial scale Qo for u and d valence quarks, as obtained in our best fits with the TMD (left 
panel) and the DGLAP (right panel) evolution, Tabic [Tlj Notice that for this analysis we have chosen to separate 
valence from sea quark contributions, while in Ref. [5] the u and d flavors included all contributions. 

This result deserves some comments. The comparison shows that the extracted u and d valence contributions, 
at the initial scale Qo — I GeV, are definitely larger for the TMD evolution fit. This reflects the TMD evolution 
property, according to which the Sivers functions are strongly suppressed with increasing Q 2 , which is not the 



12 




FIG. 5: The results obtained from the TMD-evolution fit (left panel) and from the DGLAP-evolution fit (right panel) 
of the SIDIS J 4yy'-* > ' l_ '* s ' Sivers asymmetries (red, solid lines) are compared with the HERMES experimental data [TT] 
for charged and neutral pion production. The shaded area corresponds to the statistical uncertainty of the parameters, 
see Appendix A of Ref. [5] for further details. 




FIG. 6: The results obtained from the TMD-evolution fit (left panel) and from the DGLAP-evolution fit (right panel) of 
the SIDIS A a | ^S' l ' h ~' l ' s " , Sivers asymmetries (red, solid lines) are compared with the COMPASS-p experimental data [T2] 
for charged hadron production. The shaded area corresponds to the statistical uncertainty of the parameters, see 
Appendix A of Ref. [5] for further details. 



case for the almost static collinear DGLAP evolution. Thus, in order to fit the same data at Q 2 bins ranging from 
1.3 to 20.5 GeV 2 , the TMD evolving Sivers functions must start from higher values at Qo = 1 GeV. The Sivers 
distributions previously extracted, with the DGLAP evolution, in Refs. [SJ, [T5] were given at Q 2 — 2.4 GeV 2 ; 
one should notice that if we TMD evolve the Sivers distributions on the left side of Fig. [7] up to Q 2 = 2.4 GeV 2 
we would obtain a result very close to that of Refs. [T3] (and to that of the right side of Fig. [jj) . 



III. CONCLUSIONS AND FURTHER REMARKS 



We have addressed the issue of testing whether or not the recently proposed Q 2 evolution of the TMDs 
(TMD-evolution) can already be observed in the available SIDIS data on the Sivers asymmetry. It is a first 
crucial step towards the implementation, based on the TMD-evolution equations of Refs. [7HS], of a consistent 
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SIVERS FUNCTION -TMD 



SIVERS FUNCTION - DGLAP 






x 




FIG. 7: The first moment of the valence u and d Sivers functions, evaluated at Q = Qo, obtained from our best fits of the 
^sinOh-^s) azimutlial moments as measured by HERMES p] and COMPASS [H[23] Collaborations. The extraction 
of the Sivers functions on the left side takes into account the TMD-evolution (left column of Table [TTJ) , while for those 
on the right side it does not (right column of Table Ell. The shaded area corresponds to the statistical uncertainty of the 
parameters, see Appendix A of Ref. [o for further details. 



QCD framework in which to study the TMDs and their full Q 2 dependence. That would put the study of TMDs 
- and the related reconstruction of the 3-dimensional parton momentum structure of the nucleons - on a firm 
basis, comparable to that used for the integrated PDFs. 

Previous extractions of the Sivers functions from SIDIS data included some simplified treatment of the Q 2 
evolution, which essentially amounted to consider the evolution of the collinear and factorized part of the 
distribution and fragmentation functions (DGLAP-evolution). It induced modest effects, because of the slow 
Q 2 evolution and of the limited Q 2 range spanned by the available data. The situation has recently much 
progressed, for two reasons: the new TMD-evolution [§] shows a strong variation with Q 2 of the functional 
form of the unpolarized and Sivers TMDs, as functions of the intrinsic momentum k±; in addition, some new 
COMPASS results give access to Sivers asymmetries at larger Q 2 values. 

It appears then possible to test the new TMD-evolution. In order to do so one has to implement the full 
machinery of the TMD-evolution equations in a viable phenomenological scheme. We have done so following 
Ref. [9] and the simplified version of the TMD-evolution given in Eqs. ((sj)- (j^J) - We have used them in our 
previous procedure adopted for the extraction of the Sivers functions [51 ll3lfl8] . with the same input parameters; 
moreover, we have considered also the updated HERMES [TT] and the new COMPASS [12J data. 

A definite statement resulting from our analysis is that the best fit of all SIDIS data on the Sivers asymmetry 
using TMD-evolution, when compared with the same analysis performed with the simplified DGLAP-evolution, 
exhibits a smaller value of the total x 2 > &s shown in Table m Not only, but when analyzing the partial contribu- 
tions to the total \ 2 value of the single subsets of data, one realizes that such a smaller value mostly originates 
from the large Q 2 COMPASS data, which are greatly affected by the TMD evolution. We consider this as an 
indication in favor of the TMD evolution. 

A more comprehensive study of the TMD evolution and its phenomenological implications is now necessary. 
Both the general scheme and its application to physical processes need improvements. The recovery of the usual 
collinear DGLAP evolution equations, after integration of the TMD evolution results over the intrinsic momenta, 
has to be understood. Con sider, as an example, the simple expression of the evolution of the unpolarized TMD 
PDF, as given in Eq. (27). Such an evolution describes how the TMD dependence on k± changes with Q 2 , 
but does not induce any change in the x dependence, which, at this order, remains fixed and factorized. The 
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TABLE II: Best values of the free parameters, Eq. (41 1, for the Sivers functions of u and d valence quarks, as obtained 
from our TMD-fit, TMD-analytical-fit and DGLAP-fit, at Qo — 1 GeV. The errors reported in this table correspond 
to the maximum and minimum values of each parameter in a restricted parameter space constrained by the condition 
A\ 2 = 20, corresponding to 95.45% confidence level. They correspond to the shaded area in Fig. [j] 



TMD Evolution (exact) 



TMD Evolution (analytical) 



DGLAP Evolution 



]Yf — f) 77+O.23 


N Uv = 0.75t°; 2 f 


N Uv = 0.45±g;?? 


N dv = -1.001°;™ 


N dv = -1.00±°-.«S 


iV d „ = -1.00tg;gl 


a Uv = 0.68±°;| 


a u „ = 0.821°;!* 


a Vv = 1.081°;™ 


«*, = i-ii±S:S? 


a d „ = 1-36±1;SS 


a d „ = 1.7lJ'if 


P = 3.i±*:S 


/J =4.0±*;| 


/3 =6.9«;* 


Ml = 0.40l^ 3 GeV 2 


Ml = 0.341J -j| GeV 2 


M 2 = 0.19l°;^ GeV 2 



question whether or not one can recover the usual DGLAP evolution, which changes the x dependence, for the 
integrated PDFs arises naturally at this point. A naive integration of Eq. (27) on fej_, over the full integration 
range, would give f q / p (x, Q) = f q / p (x, Qo) R(Q, Qo) which is not the correct PDF evolution. However, the k± 
integration should have upper limits which depend on x and Q 2 , and the full TMD evolution is more complicated 



than the simplified version used here, as explained at the beginning of Section I A 



We have made a safe phenomenological use of the TMD evolution equations; it is true that they induce a strong 
change in the k± dependence of the unpolarized and Sivers TMDs, leaving unchanged the x dependent shape, 
thus neglecting the collinear DGLAP evolution, but this should not be a problem. Infact, as we have shown 
explicitly in Fig. [2] (dashed curve), the collinear DGLAP evolution is negligible in the Q 2 region considered; 
Fig. |] is dr awn for x — 0.1, but a similar conclusion holds for all x values involved in the SIDIS data used in the 
paper. Moreover, the extra factors R(Q, Qo) arising in the TMD evolution, cancel out, as already explained, in 
the expression of the Sivers asymmetries. 

A fresh analysis of TMD dependent data, both in polarized and unpolarized, SIDIS and Drell-Yan processes, 
has to be carefully performed including TMD evolution from the beginning in an unbiased way. Most importantly 
so, should predictions for future high energy experiments, like the planned EIC/ENC colliders, be considered 
or re-considered. 
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